Rotational modes in molecular magnets with antiferromagnetic 

Heisenberg exchange 

O , J. Schnack 

^vq ! Universitdt Osnahriick, Fachbereich Physik 

Q^l Barbarastr. 7, 49069 Osnabriick, Germany 

(D ' 



(N 



C/3 



C^ 



CJ 



O 



X 



Marshall Luban 

Ames Laboratory & Department of Physics and Astronomy, Iowa State University 

Ames, Iowa 50011, USA 



Abstract 

In an effort to understand the low temperature behavior of recently synthe- 
sized molecular magnets we present numerical evidence for the existence of a 
rotational band in systems of quantum spins interacting with nearest-neighbor 



Q ' antiferromagnetic Heisenberg exchange. While this result has previously been 

noted for ring arrays with an even number of spin sites, we find that it also 
applies for rings with an odd number of sites as well as for all of the polytope 
►^ , configurations we have investigated (tetrahedron, cube, octahedron, icosahe- 

CD ' dron, triangular prism, and axially truncated icosahedron). It is demonstrated 

1^ ' how the rotational band levels can in many cases be accurately predicted us- 

Q^ . ing the underlying sublattice structure of the spin array. We illustrate how 

^^ , the characteristics of the rotational band can provide valuable estimates for 

Q I the low temperature magnetic susceptibility. 

PACS: 75.10.Jm, 75.40.Cx 

S 

^ . I. INTRODUCTION AND SUMMARY 

^' 

The subject of molecular magnets has greatly advanced in recent years due to notable 
progress in synthesizing bulk samples of identical molecular-size unitsiTa, each containing a 
relatively small number of paramagnetic ions ( "spins" ) that mutually interact via Heisenberg 
exchange. An especially attractive feature of many of these systems is that the intermolecular 
magnetic interactions are utterly negligible as compared to the intramolecular interactions. 

Already at this early stage it is clear that even Heisenberg systems of relatively modest 
size pose a major theoretical challenge. A stunning example is provided by the recently syn- 
thesized molecular magnetH {Mo72Fe3o}, where the 30 Fe^"*" ions (spins 5/2) occupy the sites 
of an icosidodecahedron. The total dimension of the Hilbert space for this spin system is a 
staggering 6^°, namely of order Avogadro's number, utterly precluding the calculation of the 
energy eigenvalues and eigenvectors on any imagined configuration of immense, ultra-fast 
computers. This is the context for our exploration in this article of a generic feature of the 



low-lying excitation energies of a finite number of spins interacting via antiferromagnetic 
Heisenberg exchange. With the knowledge of the low-lying excitations one can establish the 
very low temperature properties and in some cases even arrive at an estimate of the temper- 
ature range for the manifestation of essentially quantum behavior. This is illustrated in the 
present article for {Mo72Fe3o}; we arrive at an expression for the temperature dependence of 
the weak-field susceptibility at very low temperatures. We also provide an estimate for the 
temperature above which this system can reliably be described by the classical Heisenberg 
model, a far more practical theoretical platform than the corresponding quantum model. 
Indeed the classical Heisenberg model is currently being exploited so as to provide detailed 
quantitative DredictionsElQ for {Mo72Fe3o} that are being compared to the results of ongoing 
experiments.^ 

In this article we focus on a generic feature of highly symmetric geometric arrays (ring 
structures or polytopes) housing a finite number of spins interacting via antiferromagnetic 
Heisenberg exchange. Whatever the intrinsic spin of the individual paramagnetic ions, or 
the specific geometrical symmetries of the spin array, there always exist what we refer to 
as rotational modes. Besides their intrinsic interest, we show that the knowledge of these 
modes can in many cases be used to obtain good estimates of physical observables, such 
as the magnetic susceptibility at very low temperatures. We consider in the following a 
finite number, N, of quantum spins, each of intrinsic spin s, that in most examples interact 
via nearest-neighbor isotropic Heisenberg exchange. We assume that all nearest-neighbor 
pairs of spins interact with the same coupling constant, one which favors antiferromagnetic 
ordering. If S denotes the total spin operator, the operator S"^ commutes with the Hamilton 
operator, and thus we can structure the set of all energy eigenvalues according to the total 
spin quantum number S, extending up to S'max = Ns . A key point of the present work 
relates to the subset of minimal energies for the allowed values of S. We may summarize 
our findings as follows: Whatever the details of the system, this subset of minimal energies 
appears to define what we shall refer to as a "rotational band", i.e., is well approximated by 
a dependence on S of the form S{S + 1). We choose the term "rotational band" to indicate 
that this portion of the spectrum is similar to that of a rigid rotor. Similar behavior is 
commonly found also in nuclear and atomic physics. 

The occurrence of a rotational band has been noted on several occasions for an even 
number of spins defining a ring structure. The minimal energies have been describedEiO as 
"following the Lande interval rule" . However, we find that the same property also occurs for 
rings with an odd number of spins as well as for the various polytope configurations we have 
investigated, in particular for quantum spins positioned on the vertices of a tetrahedron, 
cube, octahedron, icosahedron, triangular prism, and an axially truncated icosahedron. Ro- 
tational modes have also been found in the context of finite triangular lattices of spin-1/2 
Heisenberg antiferromagnetslBlla. 

Using only the sublattice structure of the various spin arrays, which is provided by 
symmetry arguments, we are able to approximate the coefficient of the S{S + 1) dependence 
to good accuracy. Our method describes in general how this approximate coefficient can 
be deduced. We can therefore obtain an estimate of the ground state energy as well as the 
low-lying rotational excitations. It is clear that at low temperatures these minimal energies 
provide the major contribution to thermal averages. This enables us to discuss the low- 
temperature behavior of quantities such as the magnetic susceptibility without knowledge 



of the complete eigenvalue spectrum. We illustrate these considerations for the special case 
of {Mo72Fe3o}. 

The layout of this article is as follows. In Sec. |I| we present our numerical findings 
for various Heisenberg spin systems and motivate in Sec. |T| how the rotational band is 
connected to the topology of spin sites. Finally in Sec. |V| we discuss some implications of 
the rotational band on physical observables. 

II. ROTATIONAL BANDS 

The Hamilton operator for the isotropic Heisenberg model in the absence of an external 
magnetic field reads 

H = -2J^ l{u)-l{v) , \fu:s{u) = s, (1) 

{u,v) 

where J is the exchange interaction with units of energy, and J < results in antiferromag- 
netic coupling. The vector operators |(m), underlined with a tilde, are the single-particle 
spin operators with eigenvalue equations 

!(«) j I s(u) m{u) ) = s{s + 1) I s{u) m{u) ) (2) 

s^ziu) I s{u) m{u) ) = m{u) \ s{u) m{u) ) . 

The sum in (|1|) runs over all distinct interacting pairs (m, v) of spins at positions u and v. 
For a closed ring with nearest-neighbor interaction the index v would simply equal m + 1 
and the sum is understood to fulfill the cyclic boundary condition. 

A. Heisenberg square 

One of the few systems that possesses a rigorous parabolic rotational band is the Heisen- 
berg square, i.e., a ring with iV = 4. Because the Hamilton operator can be rewritten 
as 

H = -J (I' - |?3 - §24) , |i3 = 1(1) + 1(3) , §24 = 1(2) + 1(4) , (3) 

with all spin operators S^, 5*^3 and 6^34 commuting with each other and with H, one can 
directly obtain the complete set of eigenenergies and these are characterized by the quantum 
numbers S", S'13 and 5*24. In particular, the lowest energy for a given total spin quantum 
number S occurs for the choice S'13 = 'S'24 = 2s 

Es^^in = -J[S{S+l)-2-2s {2s + 1)] = Eo - J S (S + 1) , (4) 

where Eq = 4s(2s + 1) J is the exact ground state energy. The various energies Es^min form a 
rigorous parabolic rotational band of excitation energies. Therefore, these energies coincide 
with a parabolic fit (crosses connected by the dashed line on the l.h.s. of Fig. |ip passing 
through the antiferromagnetic ground state energy and the highest energy level, i.e., the 
ground state energy of the corresponding ferromagnetically coupled system. 



B. Heisenberg rings with A'^ > 4 

We have calculated all energy levels by diagonalizingpl the Hamilton matrix for a variety 
of rings with different values of N and s. All of these systems exhibit a rotational band, 
irrespective whether N is even or odd and for both integer and half-integer values of s. 
That is, the subset of minimal energies is well approximated by a dependence on S which is 
proportional to S{S + 1), i.e., follows the Lande interval ruleETM, 

Es,rmn ^ E, - J ^^^^ S{S + 1) . (5) 

We determine the parameters Ea and D{N, s) so that formula (||) reproduces our calculated 
values of the lowest and highest energies of the rotational band, i.e., the ground state energy 
of the antiferromagnetic system and the ground state energy of the corresponding ferromag- 
netic system. In all cases we have observed, that if deviations occur the fitting parabola of 
(I) lies below the rotational band. Use of (H) is illustrated in Fig. |I] (r.h.s.) for the case 
A^ = 6 and s = 3/2. The figure shows the complete spectrum (horizontal dashes) as well as 
the fit according to (|5D (crosses connected by a dashed line). One observes that the fit very 
nearly matches the energies of the rotational band, meaning that the Lande interval rule is 
obeyed with high accuracy. As a second example, spectra are shown for rings of five spins 
with s = 2 (Fig. D, l.h.s.) and s = 5/2 (Fig. ^, r.h.s.). Inspecting the low- lying excitations 
one notices that the rotational band for odd rings is not separated from the remaining states 
as much as it is for even rings. This remark also pertains to other, larger odd values of A^. 

In Table | we collect the coefficients D{N,s) for the rings we have investigated. In all 
cases D{N, s) ~ 4. For the odd rings the values of D{N, s may be somewhat smaller than 
four. We will dwell on this fact in Sec. |T|. 

It should be noted that in the large- A^-limit the rotational levels (|^) become degenerate 
since D{N,s) remains finite. Therefore, excitations within the rotational band should not 
be confused with magnons. 

C. Frustrated spin rings 

Even for spin rings with next-nearest neighbor interaction the rotational band persists. 
The energy spectra of Fig. |^ have been calculated for the Hamilton operator 

N N 

H=-2 Jnn Y^ 1{U) ■ |(m + 1) - 2 Jnnn ^^ K") " K" + 2) , (6) 

u=l u=l 

where all spins have been taken to be s{u) = 3/2. In Fig. |^ we display spectra for a ring 
with A^ = 6 and s = 3/2 for various ratios of the two coupling constants J„„ and Jnnn- 
Although details of the spectra differ, the overall appearance persists, and in particular the 
minimal energies define a rotational band. 



D. Heisenberg polytopes 

In order to illustrate the generality of the rotational band we provide several additional 
examples, the tetrahedron, the cube (Fig. |[), the octahedron (Fig. ^, l.h.s.), and the icosahe- 
dron (Fig. ||, r.h.s.). As in the previous cases the displayed energy eigenvalues are calculated 
by numerical diagonalization, except for the tetrahedron and octahedron which can be solved 
analytically. 

The tetrahedron of spins is a worthy textbook problem; it can be solved with a few lines 
of algebra because the Hamilton operator simplifies to 

H = -J {S^-Af) . (7) 

Therefore, the spectrum of this system consists exclusively of a rotational band. 

The case of the octahedron is similar to the Heisenberg square; the Hamilton operator 
can be written as 

H = -J{S'-S\-SI-SI) , (8) 

where Sa, Sb, Sc are the sums for pairs of spins situated at opposite vertices of the octa- 
hedron, and S is the total spin. The spin operators S"^, S% S% and S^ commute with each 
other and with H. Thus the eigenvalues of H may be written down at once and they are 
given in terms of the quantum numbers S, Sa, Sb and Sc- Therefore, the lowest energy for 
a given value of S is achieved if Sa = Sb = Sc = ^s, and its value is given by 

Es,rmn = - J [S (S + I) - 3 ■ 2s (2s + 1)] . (9) 

This is another one of the few cases where the minimal energies define a rigorous rotational 
band (see Fig. |], l.h.s.). 

The remaining examples of the cube (Fig. ^) and the icosahedron (Fig. ^ r.h.s.) illustrate 
the behavior of the gaps between the rotational band and the remaining eigenenergies. It is 
worth noting that the rotational band of the icosahedron is not as well separated from higher 
energy levels as it is for the cube. This behaviour is similar to that discussed above for even 
and odd rings. Systems which are bipartite, i.e., can be subdivided into two sublattices 
with interactions only between spins of different sublattices (rings with even N and the 
cube), show a significant gap, whereas systems, that are non-bipartite appear to show much 
smaller gaps." The latter systems are often also called frustrated.li^ Two other cases we 
have studied, the equilateral triangle prism and an axially truncated icosahedron, conform 
with these trends. 



III. CONJECTURE ON ROTATIONAL BANDS 

A. Heisenberg rings; even N 

It turns out that an accurate approximate formula for the coefficient D{N, s) of (^) can 
be developed using the sublattice structure of the spin array. As an introductory example 



we repeat the basic ideas for Heisenberg rings with an even number of spin sitestJ. Such 
rings are bipartite and can be decomposed into two sublattices, labeled A and B, with every 
second spin belonging to the same sublattice. From classical spin dynamics it is known 
that the classical ground state, sometimes called the classical Neel states, is given by an 
alternating sequence of opposite spin directions. On each sublattice the spins are mutually 
parallel. Therefore, a quantum trial state, where the individual spins on each sublattice 
are coupled to their maximum values, Sa = Sb = Ns/2, could be expected to provide a 
reasonable approximation to the true ground state, especially if s assumes large values. Such 
trial states are called Neel-like. For rings with even N the approximation to the respective 
minimal energies for each value of the total spin S = Sa + Sb is found to be given bylB 



4 7 

papprox ^ '-' 

S,min T\T 



S(S. l)-2^(^.l 



(10) 



This approximation exactly reproduces the energy of the highest energy eigenvalue, i.e., the 
ground state energy of the corresponding ferromagnetically coupled system {S = Ns), since 
the true eigenstate with all spins assuming their largest m quantum number, m = s, is a 
linear combination of Neel-like states. For all smaller S the approximate minimal energy 
^ strain ^^ bouudcd from below by the true one (Rayleigh-Ritz variational principle). The 
solid curve displays this behavior for the example of A^ = 6, s = 3/2 in Fig. |l] (r.h.s.). The 
entries in Table | provide additional numerical support for the approximation D[N, s) ~ 4 
adopted in ([T0|). For each fixed even N the coefficient D{N, s) approaches 4 with increasing 



s. 

The approximate spectrum, (|l3), is similar to that of two spins, Sa and Sb, each of 
spin quantum number Ns/2, that are coupled by an effective interaction of strength 4J/N. 
Therefore, one can equally well say, that the approximate rotational band considered in ( p!OD 
is associated with an effective Hamilton operator 



4 7 

rrapprox 

~ ~ N 



S' - Si - s' 



B 



(11) 



where the two sublattice spins, Sa, Sb, assume their maximal value Sa = Sb = Ns/2. 
This result agrees with that obtained bv a different procedure, which utilizes a Fourier 
representation of the Hamilton operator.li3 

In retrospect one realizes that all we needed for rings with even N to arrive at (|1^) was 
a sublattice structure in order to build Neel-like trial states. The sublattice structure can be 
deduced from the classical ground state or the symmetries of the spin array, which manifest 
themself in the classical ground state. In the case of rings with an even number of sites the 
symmetry is the cyclic shift symmetry: Sublattice A can be transformed into sublattice B 
by a single shift. These considerations lead us in the following subsection to a generalization 
for other systems. 

B. Conjecture 

We assume that the Heisenberg spin system can be decomposed into Ngi sublattices 
according to a symmetry transformation. The spins of each sublattice are to be coupled to 

6 



their maximal values Ssi = Ns/Ngi. Then we conjecture that the energies of the rotational 
band can be approximated as 

E^Zn =-J^ [S{S + 1) - iV., SsASsi + 1)] , (12) 

where the parameter D is to be fixed by the requirement that the energy -E'fcrro = —2JNb s^ 
of the corresponding ferromagnetic ground state, for which S = Ns, is reproduced exactly. 
The energy -E'fcrro depends only on the number of distinct bonds A*";,. Therefore, the coefficient 
D is given by 

which is independent of s. For rings with an even number of sites this formula reproduces 
the value D = 4, whereas for the cube, which has Nj, = 12, it gives D = 6. That these 
values of D provide a very good estimate of the coefficients D{N, s) is demonstrated in Fig. |l| 
(r.h.s.) for a ring with N = 6 and s = 3/2, and in Fig. ^ for cubes with s = 1 and s = 3/2. 
If there are several possible partitions into sublattices according to various symmetries of 
the spin array one can safely choose that one which leads to the lowest approximate ground 
state energy, because this energy is bounded from below by the true ground state energy. 
For spin arrays with low symmetry the gained approximation can be rather poor. This is 
for instance the case for rings of odd N. The classical ground state structure would suggest 
as many sublattices as spin sites, which leads to rather small coefficients D as depicted 
in Table | ("simple conjecture"). A closer inspection shows that it is possible to obtain a 
better approximation if one defines the approximate rotational band by the parabola passing 
through the classical Neel ground state energy level and through the ground state energy 
level of the corresponding ferromagnetic system. For rings of odd N the classical Neel ground 
state energy is 

E^^ei = -2 J Ns'^ cos {{N-1)tt/N), (14) 

as can be shown using spin coherent states.li3 This leads to the "refined conjecture" for 
odd rings, see Table |. For other spin arrays, such as the icosahedron, it may be simply 
impossible to derive a good approximation for the coefficient D{N,s), nevertheless, the 
minimal energies always form a rotational band according to (|^). 

Using (|12D and (0) one can find an approximation for the rotational band for larger 
polytopes such as the icosidodecahedron, which characterizes the sites of the 30 paramagnetic 
Fe^^ ions in the recently synthesized molecular magnetQ {Mo72Fe3o}. The related Hilbert 
space has a dimension of (2s + 1)^ = 6^^ which is of the order of Avogadro's number. 
Numerical diagonalization of the Hamilton matrix is totally out of reach. Nevertheless, 
we can estimate both the ground state energy and the form of the rotational band. The 
icosidodecahedron, consisting of 20 triangles and 12 pentagons, has threefold rotational 
symmetry and thus three sublattices. Assuming nearest neighbor interaction, we have Ni, = 
60. Thus, one ends up with 

E^';^:: = -^sis+i)+60Js(^s+^^ . (is) 



The threefold rotational symmetry of the spin array is also reflected by the structure of the 
classical ground state, thus it is no surprise that quantum and classical Neel ground state 
energies, 

K,rr = 60Js (s + 1/10) , Encc1 = 60Js2, (16) 

differ from each other only by 4% for the rather high s = 5/2. 

IV. LOW-TEMPERATURE THERMODYNAMIC PROPERTIES 

It is obvious that at low temperatures the rotational band energies provide the dominant 
contribution to thermal averages, especially if these energies are well separated from the 
remaining energy levels. This suggests, for example, that we approximate the partition 
function by 

Z{f])^ Y^ rf^e-^[^''-^^^(^+^)] , (17) 

where 1^5 is the degeneracy factor of the eigenenergies belonging to the rotational band, 
and Ea = J{D/N) NsiSsi{Ssi + 1) according to (|I^). For bipartite systems, i.e. systems 
that can be subdivided into two sublattices according to the theorem of the Lieb-Schultz- 
Mattisialli, ds = 25 + 1. For non-bipartite systems not much is known about the degeneracy. 
First investigations show that it is possible to establish rules for the degeneracy of certain 
states also for non-bipartite systems.c^ For the approximate rotational energies we take the 
degeneracies to be that resulting from the coupling of the sublattice spins.t3 It might be that 
the full Hamiltonian lifts this degeneracy somewhat, nevertheless this procedure provides a 
reasonable approximation for the degeneracies of the true minimal energies. 

The corresponding approximation of the zero-field susceptibility is then given by 

n Smax 1 / M^S \ 

Xo « ^ E ^ ( E mA ,-.|-.-^^^W..1 . (18) 

Inspecting Fig. |^ one sees how the rotational band contributes to the susceptibility in the 
case of the ring with N = 6 and s = 5/2. The rise at low temperatures is mostly determined 
by the first excited level of the rotational band (dashed curve labeled 1). 

As our final example we show in Fig. ^ the result using (|T^ and (|T^) for the low 
temperature behavior of the zero-field susceptibility for {Mo72Fe3o}. For this system 
ds = min{(2S'+l)^, {2S+1){76 — S)}. The susceptibility (solid curve) rises very rapidly with 
increasing temperature to the resulted for the classical Heisenberg model (dash-dot curve). 
Inspecting ([15|) one can understand that the rapid rise is due to the small energy difference 
between ground and first excited state and in particular because of the small coefficient J/5 
of the S{S + 1) term. In addition the S = 1 level of the approximate rotational band is 
9-fold degenerate and thus the rapid rise in xo commences at very low temperatures. Thus, 
for {Mo72Fe3o}, which has a nearest neighbor coupling constantQ of J/ks ~ 0.75 K, we ex- 
pect, that weak- field susceptibility measurements will confirm the rapid decrease on cooling, 
which is a genuine quantum feature and not present in the classical counterpart, only at 
temperatures below T ^ 0.15 K. 
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TABLES 
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TABLE I. Coefficients D{N, s) for various Heisenberg rings, to be used in conjunction with 
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FIG. 1. Energy spectra of antiferromagnetically coupled Heisenberg spin rings (horizontal 
dashes) . Here and in all subsequent figures the crosses connected by the dashed line always repre- 
sent the fit to the rotational band according to (0), which by definition matches both the lowest 
and the highest energies exactly. On the l.h.s the dashed line reproduces the exact rotational band, 
whereas on the r.h.s. it only approximates it, but to high accuracy. The solid line on the r.h.s. 
corresponds to the approximation of (IH 
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FIG. 2. Energy spectra of antiferromagnetically coupled Heisenberg spin rings with A^ 
s = 2 (l.h.s.) and s = 5/2 (r.h.s.). 
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FIG. 3. Energy spectrum of an antiferromagnetically next neighbor coupled Heisenberg ring 
with competing next-nearest neighbor interaction. 
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FIG. 4. Energy spectrum of antiferromagnetically coupled Heisenberg cubes. The solid lines 
correspond to the approximation of (12) and (13), i.e. D = G. 
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FIG. 5. Energy spectrum of an antiferromagnetically coupled Heisenberg octahedron (l.h.s. 
and icosahedron (r.h.s.). 
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FIG. 6. Zero- field magnetic susceptibility of the antiferromagneticahy coupled Heisenberg 
spin ring with A^ = 6 and s = 5/2. The solid curve displays the full quantum solution; the dashed 
curves correspond to approximation (^) upon including the contributions of rotational levels up 
to 5 = 1, 5 = 3, and S = 15. 
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FIG. 7. Zero- field magnetic susceptibility of {Mo72Fe3o} calculated using (|^) and (18). The 
dashed curve is obtained using only the two lowest levels and the solid curve using all levels of the 
approximate rotational band. The classical result is given by the dash-dot line. 
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